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ABSTRACT 

We present X-ray imaging and spectral analysis of all microflares the Reuven Ramaty High Energy 
Solar Spectroscopic Imager {RHESSI) observed between March 2002 and March 2007, a total of 25,705 
events. These microflares are small flares, from low GOES C Class to below A Class (background 
subtracted) and are associated with active regions. They were found by searching the 6-12 keV 
energy range during periods when the full sensitivity of RHESSI's detectors was available (see paper 
I). Each microflare is automatically analyzed at the peak time of the 6-12 keV emission: the thermal 
source size is found by forward-fitting the complex visibilities for 4-8 keV, and the spectral parameters 
(temperature, emission measure, power-law index) are found by forward fitting a thermal plus non- 
thermal model. The combination of these parameters allows us to present the first statistical analysis of 
the thermal and non-thermal energy at the peak times of microflares. On average a RHESSI microiiaTe 
has a fitted thermal loop width 8 Mm (11"), length 23 Mm (32") and volume 1 x 10^^ cm"^, temperature 
13 MK, emission measure 3 x 10"*^ cm~^ and density of 6 x 10^ cm~'^. There is no correlation between 
the loop size and the flare magnitude, either flux in the loop or GOES class, indicating that microflares 
are not necessarily spatially small. There is also no clear correlation between the thermal parameters 
except between the RHESSI and GOES emission measures, the GOES values are generally twice the 
RHESSI emission measures. The microflare thermal energy at the time of peak emission in 6-12 keV 
ranges over 10^^ to 10^° erg and has a median value of 10^^ erg. The frequency distribution of the 
thermal energy deviates from a power-law at low and high energies arising from a deficiency of events 
due to instrumental and selection effects. It is difficult to compare this energy distribution to previous 
thermal energy distributions of transient events, as the work sought nanoflares through imaging in 
EUV or soft X-rays and covered just a few hours. There are large uncertainties in the majority of 
the non-thermal parameters, due to the steep spectra down to low energies. We typically find a 
power-law index of 7 above a break energy of 9 keV, which corresponds to a low-energy cut-off in the 
electron distribution as low as 12 keV. The resulting non-thermal power estimates, covering 10^^ to 
IQ'^s gj.g g-i ^[i]^ median value of 10^^ erg s~^, therefore have large uncertainties as well. The few 
microflares with unexpectedly large non-thermal powers 10^^ erg s~^ have the smallest uncertainties, 
of about 10%. The total non-thermal energy however is still small compared to that of large flares as 
it occurs for shorter durations. 

Subject headings: Sun:flares - Sun: corona - Sun: X-rays, gamma rays - Sun: activity 



1. INTRODUCTION 

The solar corona exhibits a myriad of transient 
energy releases over many scales, from large flares 
down to nanoflares. The frequency distribution of the 
energy in these eve n ts has been st u died extensively 
(iCrosbv et all Hool IShiiiii^ [l995l : iKrucker fc Be^ 



Aschwa.nden et al.l l2000t Parnell fc Jupr 



iLin et all 120011 : lAschwanden k Parnelll l2002l : 



Benz fc Krucked l2002l ) and has been found to be well 
represented by a power-law of the form 

dA^ = AW-"dW (1) 

where dA'^ is the number of events per unit time with 
energy between W and W + dW. These distributions 
are of particular interest as they elucidate the amount 
of energy available, how often it is released and in which 
events and form (thermal or non-thermal) it predomi- 
nantly occurs. The energy release observed in normal- 



size flares is not sufficient to constantly and consistently 
heat the corona to the observed few million Kelvin. So 
the question is then whether this could be achieved by ex- 
tending the observed flare-like energy releases to smaller 
scales. This concept can be expressed in terms of the 
power-law index of this distribution: if a > 2 then the 
smallest events have a high occurrence rate and their en- 
ergy dominates over larger flares, possibl y matching the 
energy in coronal heating (jHudsonlll991h . This requires 
the assumption that these distributions be continuous 
into the unobservable low-energy range, which is diffi- 
cult to determine as there are instrumental and selection 
effect s that both cut-off and bias the observed distribu- 
tion (|Aschwanden fc Parnelil2002f ). 

The instantaneous thermal energy in these events may 
be calculated from 



Wt = Sn^kBTV 



(2) 
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where ric is the electron density, V the volume of the 
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obtained by imaging the events, however there can be an 
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overestimate in this observed volum e to the true volume 
by a filling factor / fti 1 to 10"^ (ICargil l fc KlimchuS 
119971: iTakahashi fc Watanabell2000D . In this work we as- 
sume / = 1, which will be discussed later. The tem- 
perature and emission measure may be found cither di- 
rectly from the spectrum or by imaging with different 
wavelength filters. Assuming constant density, the emis- 
sion measure is related to the density and volume as 
EM = nlV and so the thermal energy is 



Wt = 3VEM ■ VksT. 



(3) 



Note that as losses are not taken into account here, 
the energy going into the thermal plasma will be 
larger. For the smallest events this thermal en- 
ergy has been found using pixelated detectors in EUV 
and soft X-rays, with simultaneous pixel brighten- 
ings within some area being registered as an event 



20001 Parnell fc Judd 


200011 Aschwanden & Parnellll2002l: 


Benz fc Kruckeiil2002 


). The area inferred from these, of- 



ten spatially discontinuous, brightened pixels gives an 
estimate of the volume, and observations using different 
filters give temperature and emission measure informa- 
tion. 

The events observed in EUV are termed "nano" -flares 
as they have about 10~^ times the energy in large 
flares, th e limit of their observa bility being down to 
10^"* erg (lAschwandeii et al]|2000f) . Parker's hypothet- 
ical nanofiare (|Parkeil 11981 is an estimate of the ba- 



sic unit of a localized impulsive burst of energy release, 
with energies < 10^* erg, with ensembles of them con- 
stituting the observed events. The thermal energy of 
these events outside of active regions ha s been investi- 
gated in soft X-rays with Yohkoh/SXT (jKrucker et al.l 
[l99l finding energies of 10^ ^ erg per even t, and 
in EUV using SOHO/EIT (jKrucker & Benz IMS'; 
Benz fc Kruckerl [2002.) and T RACE ([ParneU fc Jupfl 



20001 : lAschwanden et al.l l2000t ) providing energies be- 



tween lO^'' to 10^^ erg. Small events in active regions, 
termed "active region transient bri ghtenings" , we re seen 
in soft X-rays with Yohkoh/SXT ijShimizul fl995h . with 
energies between 10^^ to 10^^ erg. All of these studies 
found the power law index of the frequency distributions 
to be between a = 1.5 — 2.6. lAschwanden fc Parnell 
()2002f ) investigated the effect of instrumental bias on 
the different indices from S OHO/EIT, TRACE and 
Yohkoh/SXT data. iParnelll (j2004) later pointed out 
that the discrepancy between power-law indices from dif- 
ferent instruments can be due to the effects of least- 
squares fitting some the of binned histograms. Similar in- 
dices were obtained when a maximum likelihood method 
(jParnell fc Jupp 2000) was used instead. 

The non-thermal hard X-ray emission is assumed to 
be due to a power-law distribution of electrons emitting 
hard X-rays via bremsstrahlung in a thick target. The re- 
sulting powe r-law photon spectrum reflects this electron 
distribution (jBrownl [l971l ) allowing the power in these 
accelerated electrons above a low-energy cut-off Eq (in 
keV), to be calculated as 



Pn(>£^c) = 9.5x 102V(7-1) 



2'2)^o£^c " ergs 



(4) 



where 7 and /q are the index and normalization of the 
photon power-law spectrum (in units of photon fl ux, s~^ 
cm~^ keV~^) a nd /3(to, n) is the beta function (jBrownl 
II97II: iLinI 11974) . Therefore observing the hard X-ray 
spectrum of these events for various time intervals during 
each flare is sufficient to obtain an estimate of the non- 
thermal energy. However there is ambiguity in the low- 
energy cut-off, Ec, because the observed photon spec- 
trum depends only weakly on it, with a resulting ffatten- 
ing of the photon spectrum below ee not uniquely related 
to Ec, with eb < Eq (Holman 2003). Uncertainty in Ec 
results in a larger uncertainty in the power estimate: a 
factor of 2 increase/decrease in Eq would result in a fac- 
tor of 8 decrease/increase in the power for flat spectra 
(7 = 4) or a factor of 64 decrease/increase in the power 
for steep spectra (7 = 7). 

Previous statistical studies of the non-thermal energy 
in flares used the energy threshold of the instrument as 
an estimate of Eq , as they did not observe to low enough 
energies nor had sufficient ener gy resolution t o obs erve 
the flattening of the spectrum. ICrosbv et al.l (|1993f ) us- 
ing SMM/HXRBS estimated the non-thermal energy in 
large flares > 25 keV, finding that the power distribu- 
tion at the peak time of emission had a power-law with 
a — 1.67 over 10^^ to 10"^° erg s"""^ and the total non- 
thermal energy of these events had a = 1.53 over 10^^ 
to 10'^^ erg. Microfla res were observe d down to 8 keV 
with CGRO/BATSE (|Lin et al.ll200lb . finding energies 
over 10^7 t o 10^° erg, arid hy WATCH/GRANAT down 
to 10 keV (ICrosbv et al.l[l998h . 

The Reuven R amaty High En ergy Solar Spectroscopic 
Imager RHESSI (|Lin et al.l 2002f) is uniquely suited to in- 
vestigate both the energy in the heated and accelerated 
electron populations, thermal and non-thermal emission, 
of these small events. This is due to its unprecedented 
sensitivity to 3-25 keV X-rays, hi gh spectral resolution 
and imaging capabilities (jKrucker et al.ll2002l) . Previous 
studies of RH ESSI microflares have concentrated on indi- 
vidua l events (jKrucker et al.l 120021 : iBenz fc Grigiil2002l 
l2003f ) or small samples, often compared to other wave- 
lengths, (iLiu et al. 2004: Qiu et al. 2004; Battagha et al] 
[2005l:lKundu et al.ll2005L 120061: IStoiser et al. 2007). 

Here we present the first analysis of all RHESSI mi- 
croflares found as transient bursts in 6-12 keV during 
periods of shutter- out observ ations, as detailed in part I 
of this article (Chri ste et aIH 2008). Between March 2002 
and March 2007 25,705 events were found. These are 
active-region phenomena of low C GOES class to below 
A-class. For the analysis presented here we measure the 
energy for 16 seconds around the time of peak emission in 
6-12 keV. The peak time is only used as it presents the 
best opportunity to obtain enough counts above back- 
ground to permit analysis in each event. Given the stan- 
dard flare time profile of a sharp impulsive rise followed 
by slower decay phase, this means that this analysis will 
mainly cover the impulsive phase of these microflares. 

The calculation of the thermal energy using equa- 
tion ([3|) requires knowledge of the volume of the ther- 
mal source and thermal parameters from the spectrum. 
Imaging these events, and hence estimating the volume 
of thermal emission, is detailed in fj2l In fj3]the spectral 
fltting of the events is described, obtaining the thermal 
and non-thermal parameters of these events. The rela- 
tionship between these parameters and those found using 




Fig. 1. — (Left) Visibility forward fit model loop shape (solid line) onto the 4-8 keV visibilities amplitudes (diamonds with error bars) as 
a function of the phase for each grid. The amplitudes for each grid/subcoUimator are shown between the dashed vertical lines, indicating 
the amplitudes measured as the grids rotate from a position angle of 0° to 180° . Thus an amplitude measured in subcoUimator 3 when 
rotated to a position angle of 45° would be plotted against 3.25 on the x-axis. (Right) Resulting image of the visibilities shown in left panel, 
background image found using MEM NJIT algorithm (Schmahl ct al. 2007), overplotted contours (25 %, 50 % and 75 %) representing the 
model shape. The loop length I and width ui are quoted. Color version available in electronic edition. 



GOES data are discussed in §3.31 The calculation of the 
thermal energy and power in non-thermal electrons and 
the resulting frequency distributions is presented in fj4l 
We discuss these results further, including how this anal- 
ysis at peak time relates to the emission over the whole 
of the microflare, in fj^l 

2. IMAGING USING VISIBILITIES 

RHESSI imaging is achieved through a Fourier-based 
method using r otatio n modulation collimators (RMCs) 
(jHurford et all l2002t l. Each RMC time- modulates 
sources whose size scale is smaller than its resolution. 
This spatial information encoded in the time-modulation 
profile is normally reconstructe d into an image via t ech- 
niques such as back projection (jHurford et al.ll2"002f ). 

A recently implemented alternative technique converts 
the time- modulation profile to com plex visibilities befor e 
recovering the spatial information (jHurford et al.|[2005f ). 
Each RHESSI visibility is a calibrated measurement of a 
single Fourier component of the source distribution mea- 
sured at a specific spatial frequency, energy and time 
range. The visibilities are the complex quantities ob- 
tained from the fitted amplitude and phase of the mod- 
ulated time profile for a particular roll orientation (ro- 
tational phase) of an individual RMC. The resulting set 
of visibilities for all roll angles and RMCs is a calibrated 
and compact representation of the original time profile, 
with little loss of information. 

The advantage of using visibilities is twofold. First, as 
a smaller data set has to be processed, the image recon- 
struction from the visibilities is considerably quicker than 
using the time profile directly. Second, the visibilities are 
fully calibrated measurements, representing an interme- 
diate step between the modulation profile and imaging, 
meaning that spatial information can be found directly 
from the visibilities without having to compute an image. 
This has been implemented in a Visibility Forward Fit 



(VFF) algorithm (jHurford et al.ll2005h which determines 
the best-fit parameters, with statistical errors, for simple 
assumed source geometries (elliptical Gaussian, curved 
elliptical Gaussian, multiple sources etc). 

We are primarily interested in the size of the thermal 
source, in order to make an estimate of the density and 
hence thermal energy, see ^ The images of thermal 
sources in microfiares are taken over 4-8 keV and gener- 
ally have a single elliptical source or loop shape. We fit a 
2D model of a curved elliptical Gaussian to the 4-8 keV 
visibilities for 16 sees about the peak in 6-12 keV for 
each microflare. This attempts to fit a Gaussian profile 
along the curved semi-major axis, equivalent to the loop 
arc length. If the source is not appreciably curved an 
elliptical geometry of zero curvature is returned. Seven 
parameters are obtained from each fit of this 2D model: 
the centroid position (x,y), photon flux, FWHM loop 
length and width, curvature and position angle of the 
semi-major. 

An example of this fitting is shown in Figure [1] Here 
in the left panel we have the visibility amplitudes, with 
statistical errors, for each grid and is position angle. The 
solid line is the VFF model loop shape which has fitted 
subcoUimators 3,4,5,6,8 and 9 well. For this event the 
finest subcoUimator, 1, is dominated by noise and so has 
little infiuence on the fit. SubcoUimators 2 and 7 are 
not included as they provide a poor response in this en- 
ergy range. The right panel shows the resulting image of 
these vi sibilities, produced using the MEM NJIT algo- 
rithm (jSchmahl et al.ll2007f) . with the VFF model loop 
from the left panel overplotted as contours. The contours 
correspond well with the background image. Further ex- 
amples of the resulting fits are shown in Figure [21 again 
calculated using detectors 1,3,4,5,6,8 and 9. 

The model fit for each of the microfiares not only pro- 
vides the spatial information but also a measure of the 
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Fig. 2. — Example microflares showing the MEM NJIT image (background) and the forward fit model shape (foreground contours of 
25%, 50% and 75%). The model loop length I and width w are quoted in each panel. Color version available in electronic edition. 
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Fig. 3. — Histogram of the count rate per detector in 4-8 keV, 
the different lines indicating all microflares (light grey dotted) and 
subsets with a good (solid black) and bad (solid grey) model fits 
to the visibilities. 

"quality" of the fit, based upon whether the fit con- 
verged, whether the fit parameters reach the hmit of their 
range and the size of the errors relative to the parame- 
ter. Such an objective measure of the fit quality is vi- 
tally important for an automated analysis project as it 
is impractical to visually inspect over 25,000 images. Af- 
ter processing all microflares we have 18,656 microflares 
to which the model achieved a satisfactory fit and were 
resolved, i.e. returning spatial sizes larger than the in- 
strumental resolution, 2.3". The majority of events pro- 
ducing poor fits were those that had the fewest counts. 
This can be seen in Figure [3] where the histogram of the 
4-8 keV count rate per detector of all the microflares 
is shown, as well as for the subsets of events producing 
good and bad fits. There are some microflares with large 



count rates but poor fits. This is likely due to an in- 
strumental issue, such as to an absence of spacecraft roll 
information. 

The histograms of the loop FWHM arc length I and 
width w (at loop mid-point) for the events with good 
VFFs are shown in Figure U] We find that the median 
FWHM loop arc length is 31.6" (23 Mm) and width is 
10.5" (8 Mm). The lengths distribution has a sharp peak, 
symmetrical in log-space, away from the resolution limit. 
The histogram of the ratio of the loop arc length to width 
(middle panel of FigureS]) shows that the majority of the 
microflare thermal sources are elongated structures, with 
the median value of the arc length being 3 times the loop 
width. The size of these loops shows no correlation with 
the magnitude of the flare, either the flux in the loop 
or the background subtracted GOES class (Figure [5]). 
This shows that small flares are not necessarily spatially 
small, which is certainly the case for the examples shown 
in Figures [T] and O 

The volume of this thermal emission can be estimated 
by assuming that the observed 2D loop structure has a 
cylindrical geometry as 



V = TT 



(5) 



where I is the FWHM loop arc length and w is the width 
at the loop mid-point. The histogram of the loop volume 
for the good events is shown in the right panel in Figure 
m The median volume is about 1 x 10^^ cm^'^, which 
is a factor of «66 larger than the minimum measurable 
volume of 1.5 x 10^^ cm""^, found by taking w = I — 2.3" 
in equation ([5]). 

We obtain other useful parameters from VFF, in par- 
ticular a measure of the total 4-8 keV photon flux 
from the loop. Since unmodulated background does not 
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Fig. 4. — {Left) Histogram of the model loop arc length and width for 18,656 RHESSI microfiares. The dotted vertical line indicates 
the instrumental resolution limit of 2.3". (Middle) Histogram of the ratio of 4-8 keV microflare loop length to width. {Right) Histogram 
of the microflare volume of the 4-8 keV loops, assuming cylindrical geometry from the 2D fitted width and length. The dotted vertical line 
indicates the minimum measurable model volume. 
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Fig. 5. — Correlation plots of the model loop arc length against 
the background subtracted GOES 1-8 A soft X-ray flux {top) and 
the RHESSI image model's 4-8 keV flux {bottom). 

affect the visibilities this flux measure is intrinsically 
background-subtracted as it is the emission from only 
the loop. In Figure [6] we have the differential frequency 
distribution of this 4-8 keV photon flux. This distribu- 
tion covers a range of 5 to 5000 photons s~^ cm~^ kcV^^ 
and power-law para meters were found usin g a maximum 
likelihood method (jParnefl fc JuppI l2000f ). This tech- 
nique uses the standard statistical procedure of the max- 
imum likelihood estimation to flt a skew-Laplace distri- 
bution to the data. This distribution consists of a bro- 
ken power-law; the index above the break is the true 
distribution, whereas below it the power-law fits the 



Fig. 6. — Frequency distribution of the total photon flux mea- 
sured from the VFF 4-8 keV loops. The power-law fit is performed 
using a maximum likelihood method of Parnell fc Jupp, (.2000.) and 
is not a graphical fit to the histogram shown. 

flattening/turning-over of the distribution from under- 
sampling the smallest events. So from a simple calcula- 
tion on the sample (in this case 4-8 keV fluxes) an ob- 
jective measure is obtained of the power-law index above 
a break (with errors found from the 95% sample confi- 
dence) instead of subjectively choosing bin sizes before 
line fitting a histogram. For simplicity the resulting fit 
is shown overplotted to a standard histogram in Figure 
[6l Over two orders of magnitude the power-law index is 
1.71 ± 0.02. This is steeper than the index of 1 .59 found 
for hard X-rays > 25 keV (|Crosbv et al.lll993f) b ut flat- 
ter than the that found for soft X-ra ys, 1. 7 - 2.1 (iDrakel 



19711: iLee et"allll995l:ll^ldman et al.](l997HVeronig et al 



2002f ). The distribution in Figure [6] deviates from a 
power-law at both low and high fluxes due to instru- 
mental selection effects. The events with the smallest 
and largest fluxes are missing as we are unable to suc- 
cessfully analyze these events: the smallest are hard to 
observe above background and the largest have excessive 
counts causing high detector deadtime or are excluded 
from our microflare list as RHESSI' s attenuating shut- 
ters were deployed. 

Another way in which we can use this RHESSI 4^8 keV 
thermal flux is by comparing it to the emission observed 
by GOES in its 1-8 A band (Figure [T]). Here we see 
that there is a power-law correlation, with index close 
to one, although there is some spread about the fitted 
fine. This suggests that RHESSI and GOES are observ- 
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Fig. 7. — Background-subtracted GOES 1-8 A soft X-ray flux 
against the RHESSI image model's 4-8 keV flux. 
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Fig. 8. — Histogram of the ratio of flare signal to pre-flare back- 
ground in 6-12 keV (dashed black), 4-8 keV (solid light grey) and 
12-16 keV (dashed-dotted grey) at the time of peak emission in 
6-12 keV. 

ing the emission from the same thermal plasma with the 
different temperature distribution from event to event 
accounting for the spread in the correlatio n. In compar- 
ison to RHESSI quiet Sun observations l| Hannah et al.l 
|2007( ). the smallest microflare flux measured here is over 

2 orders of magnitude larger than the limit found from 
an active-region-free quiet Sun. The quiet Sun RHESSI 
flux also correlates with the GOES flux with a slightly 
steeper slope (1.08 ± 0.13) than the microflares shown 
here. 

3. SPECTRUM FITTING 

The spectrum of each microflare is determined over 
the same time period as the visibilities in Sj2l 16 seconds 
around the time of peak emission in 6-12 keV. These 
spectra are made with 1/3 keV energy bins between 

3 keV to 30 keV using detectors 1,3,4,5,8 and 9. For 
each microflare, a background interval before and after 
the event were determined automatically with the selec- 
tion of the se background tim es described in Part I of 
this paper (jChriste et al.ll2008D . Only the pre-flare back- 
ground is used in the subtraction from the spectrum, as 
it is sharply defined by the microflare's impulsive phase. 
This was possible for 19,441 microflares. The histogram 
of the flare signal to background ratio is shown in Fig- 
ure [3 for the 6-12 keV, 4-8 keV and 12-16 keV energy 



bands. Requiring the flare signal to be at least 3 times 
the pre-flare background, we find 12,814 events suitable 
in 6-12 keV (the energy range the events were found in) 
12,472 in 4-8 keV from predominantly thermal emission 
and 9,681 in 12-16 keV from mostly non-thermal emis- 
sion. This indicates that the thermal component for more 
microflares can be obtained than the non-thermal com- 
ponent. This does not mean that many microflares do 
not have a non-thermal component, just that we cannot 
distinguish it from the background in these cases. 

These background-subtracted observed count spectra 
are forward-fitted with a model in photon space con- 
verted back to count space using the ful l RHESSI detec- 
tor response matrix (jSmith et al.l I2002D in the OSPEX 
software p ackage, an updated version of the SPEX code 
(jSchwartd [1996) . The model has both a thermal and 
non-thermal component so that we can recover the re- 
spective parameters to calculate the energy in both the 
heated and accelerated electrons. The thermal compo- 
nent contains the isothermal bremsstrahlung emission 
(free-free and free-bound) as well as hne emission from 
the C HIANTI database (|Dere et al.l 119971: iLandi et al] 
[200^ . This emission depends on the temperature T and 
emission measure EM = rilV . The non-thermal compo- 
nent is assumed to be thick-target emission of a power- 
law distribution of electrons above a low energy cut-off 
Eq, with th e photon spect rum found through numerical 
integration (|Holmanll2003D . Although this numerical in- 
tegration provides an accurate representation of the non- 
thermal emission it is too slow to compute in this fitting 
procedure, requiring multiple iterations per microflare, 
with tens of thousands of microflares to fit. Instead the 
non-thermal component is fitted with a broken power- 
law, which has an index of —7 above the break energy of 
ee and a fixed index of —1.5 below the break. An exam- 
ple of this approximation to the numerical integration is 
shown in the left panel of Figure [20] and the relationship 
between these two models is detailed further in > j4.2[ as it 
is important for calculating the power in the non-thermal 
electrons via equation (jl]). 

The fit to each microfiare spectrum is conducted us- 
ing the following strategy. First the thermal parameters 
are varied to fit the spectrum over 4-8 keV, where the 
thermal emission normally dominates. Then the fit is 
repeated but the non-thermal parameters are allowed to 
vary, while keeping the previously found thermal param- 
eters fixed, to fit the spectrum from 10 keV up to ei- 
ther 30 keV or to where the background dominates over 
the flare signal. The fit is repeated a flnal time allow- 
ing all the fit parameters to vary to fit the spectrum 
from 3 keV up to either 30 keV or to where the back- 
ground dominates over the flare signal. Examples of typ- 
ical microflare spectra and the flts are shown in Figure [9l 
The microflares here illustrate similar characteristics as 
seen in previous RHE SSI microflare studies, for example 
iKrucker et all (|2002[ ). At low energies (< 10 keV) the 
thermal component dominates with the expected spec- 
tral lines, Fe K-shell feature (about 6.7 keV) and the 
Fe/Ni lines (a bout 8 keV), for this temperature range 
(|Phillipsll200l) . At higher energies (> 10 keV) there is 
a power-law component that dominates over the ther- 
mal model which is normally assumed to be the non- 
thermal emission. Although this could be an additional 
hotter thermal component several other arguments im- 
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Fig. 9. — Example microflare spectra. Shown are the background subtracted data (black crosses), thermal model (solid line, red in 
electronic edition), non-thermal model (dashed line, blue in electronic edition) and total model (solid line, green in electronic edition). The 
title shows the GOES class and the background-subtracted equivalent. 



ly non-thermality: the presence of the Neupert effect 
in some case s, imaged hard X- 
ray footpoints (|Krucker et al.1 12002D . a nd complemen- 
tary radio and microwave observations (iLiu et al.1 120041 : 
iQiu et al.l[200l iKundu et al.1 12005'. '2006)! 

As with the imaging, the spectral fitting returns pa- 
rameters that objectively measure the quality of the 
fit, such as the fit oi' whether the fit parameters 
have large errors or reach the limit of their chosen fit- 
ting range. We obtain 9,161 events for which we trust 
the fit to the thermal component of the spectrum and 
the forward fit model loop. This is out of a possi- 
ble 9,693 microflares with good background subtraction, 
flare signal-to-background over 4-8 keV and a successful 
model fit. For the thermal and non-thermal fit as well as 
the VFF model loop fit to the visibilities we obtain 4,236 
trustworthy events. This is out of a possible 8,046 mi- 
croflares with good background subtraction, flare signal- 
to-background over 10-12 keV and successful model fit. 
The histograms of these fitted parameters are shown in 
Figure [TO] and discussed in ^ 33. II and ^ 33. 21 

3.1. Thermal Parameters 

The histograms of the fitted temperature and emission 
measure for the events with good background subtrac- 
tions and thermal fits, 9,161 microfiares, are shown in 
the top row of Figure [TOl This is about a third of the 
total sample but shows nearly all of the events with good 
background subtraction, flare signal-to-background and 
visibility forward fits (9,693 microflares). The majority of 
the temperatures found lie within a tight range of about 
10 to 15 MK, with the median temperature of around 
13MK. The emission measures vary considerably more 
than the temperatures, having a range covering over two 



orders of magnitude between 10 to 10 cni^-^. The 
median emission measure is 3 x lO''^ cm"'^. Also shown 
in Figure [10] are the average ratio of the error in the 
fit to the fitted parameter. For the temperatures this 
statistical error in the fit is < 1% and is approximately 
constant for the temperatures found. The error in the 
emission measure is w 10% at 10^^ cm"-^ but drops to 
« 1% at 10**^ cm^'^. The larger relative error in the 
events with smallest emission measure is due to the emis- 
sion measure being directly proportional to the thermal 
emission model, and so are events with small noisy spec- 
trum. The range of these parameters is discussed further 
in > j3.3l With these emission measures and the volumes 
of the emitting plasma (see ^J2]) an estimate of the elec- 
tron density can be made. The histogram of these den- 
sities is shown in Figure [TT] and range from 6 x 10^ to 
6 X lO"'^" cm^'^ with median value of 6 x 10^ cm~'^. This 
is larger than typic al coronal conditions but reasonable 
for a flaring loop (jPhilhps et al.l Il996t iGallagher et al.1 



With these fitted thermal components we can estimate 
the 4-8 keV flux from these spectrum fits and compare it 
to the flux derived from the imaging in fj^l This provides 
a consistency check to verify that these two vastly differ- 
ent analysis techniques recover similar fluxes. The his- 
togram of the ratio of the flux found from the spectrum 
model to the image value is shown in Figure [T^ The me- 
dian of these is 1.0 with some spread about this value. 
This is expected as different detectors were used for imag- 
ing and spectral analysis (additional use of detector 5 in 
imaging) and the imaging calculation uses only the di- 
agonal elements of RHESSI's detector response matrix, 
(jSmith et al.l l2002f ) , whereas the spectrum fitting uses 
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Fig. 10. — (Top) Histogram of the microflare temperature {left) and emission measure (right) from the spectral fits of 9,161 RHESSI 
microflares. {Bottom) Histogram of the photon power-law index 7 {left) and break energy {fight) for 4,236 RHESSI microflares. The 
grey lines in each panel refer to the righthand axis and show the average ratio of the error in the fit to the fitted parameter, as a function 
of the parameter. 
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Fig. 11. — Histograms of the density of the 4-8 keV loops, found 
from the emission measure and thermal volume for 9,161 RHESSI 
microflares. 

the fuU response matrix. 

3.2. Non-thermal Parameters 

Histograms of the index 7 and break energy ee of the 
non-thermal broken power-law are shown in the bottom 
row of Figure 1101 The power-law index 7 has values 
mostly ranging over 4 to 10 with the median about 7. 
This is considerably steeper than large flares observed 
Lis sed in pr 

DUqqI; H 

break energy eb ranges over 7 keV to 12 keV with the me- 
dian being a bout 9 keV, which is smaller than is found for 
larger flares (jSaint-Hilaire fc Bendl200"5h . In larger flares, 
the lower energy non-thermal emission would be masked 
by the thermal emission to tens of keV. The steep power 
laws starting at low energies leads to a strong selection 



by RHES SI, as discus sed in previous RHESSI microflare 
work f)KriIcker et all I200I : iBenz fc Grigis 2002. ). The 
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Ratio Spectrum/Image, 4-8 l<eV 

Fig. 12. — Histogram of the ratio of the 4-8 keV fiux found from 
the fitted spectrum model and the image model loop. 

effect. This is because such steep power-laws, e xtending 
down to energies where there are spectral lines (Phillip^ 
[200l . are difflcult to distinguish from a thermal compo- 
nent. A conservative approach has been taken here to 
remove any events where there is an ambiguity between 
the thermal and non-thermal components, so we discount 
any events any events with ee < 7 keV. The result is that 
we have only 4,236 microflares, about a fifth of the total 
sample. 

Also shown in Figure [10] are the average ratios of the 
error in the fit to the fitted parameter. For the power-law 
index 7 this statistical error in the fit is « 4% for 7 < 8 
and increases to « 10% for 7 > 8. This increase in the 
error shows the greater uncertainty in trying to fit steep 
spectrum. The error in the break energy is « 10% at 
about ee — 7 keV and decreases to > 1% by ee = 12 keV. 
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Fig. 13. — Microflare temperature against emissio n measure. The corre lation plot is the microflares in this study with black crosses 
representing the result s of a previous RHESSI study I IBattaelia et al.ll2005l ). The dashed straight line represents the correlation found by 
IFeldman etaTI Ill996bl) from GOES and Yohkoh/BCS. The numbered curved lines are of constant count rate per detector in 4-8 keV for 
the thermal model, as a function of temperature and emission measure. The white contour show the events that have a GOSS temperature 
< 10 MK in Figure [Til 

The large errors at low break energies shows the greater 
uncertainty in trying to separate the thermal and non- 
thermal components below 10 keV. These non-thermal 
parameters are discussed further in 
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Fig. 14. — Correlation plots of temperature (top) and emission 
measure (bottom) from RHESSI and GOES for the same peak 
time in RHESSI 6-12 keV for 6,740 microflares. The dashed lines 
indicate a one-to-one correlation. 



3.3. Correlation Between Parameters 

Figure [13] plots the temperature against emission mea- 
sure for the RHESSI microflares. There is no clear cor- 
relation between these parameters. Also shown are num- 
bered contours indicating constant count rate per de- 
tector over 4-8 keV for the thermal model as a func- 
tion of temperature and emission measure. All of the 
microflares lie between the 10 and 10'* counts con- 
tours, consistent with the non-background subtracted 
count rates for good fits shown in Figure [31 Although 
any temperature and emission measure between these 
contours could be expected, the temperatures lie in a 
tight range, mostly between 10 MK and 15 MK, with al- 
most all possible emission measures, from 10"*^ cm~'^ to 
10*^ cm~^, for this temperature range found. The model 
of the thermal emission is directly proportional to EM 
and increases with larger T, though not directly, with 
the continuum rising and flattening and the line features 
becom ing more prominent (Tandbcrg-Hansscn & Emsli3 
11988( 1. This results in the errors in the temperature and 
emission measure being anti-correlated. The thermal 
model also includes spectral features and they provide 
additional emission, particularly from the Fe K-shell over 
4-8 keV, for temperatures above 8 MK (PhiUips 2004). 
The fact that only temperatures above 8 MK have been 
found is more suggestive of a selection effect primarily 
affecting the temperatures and not a physically signifi- 
cant discovery of microflares with a lack of low tempera- 
tures and high emission measures. This selection effect is 
consistent with RHESSVs sensitivity being temperature- 
dependent. The combination of this greater sensitivity 
to hotter plasma and the differential emission measure 
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Fig. 15. — Microflare temperature (top) and emission mea- 
sure (bottom) against the background subtracted GOES l-SA fiux. 
The correlation plots are the microflares in this study, with the 
dot-dashed line (bottom panel) the fit to this data. The black 
crosses and dashed line were found from a previous RHESSI study 
jBattaglia ct al. 200 5). The solid l i ne (top panel) represents the 
correlation found bv lFeldman et aLl Ill996bl ). 

DEM decreasing as the temperature increases could ex- 
plain the tight range of temperatures found with the peak 
of this sensitivity for RHESSI shutter-out mode occur- 
ring in this temperature range. Note that this selection 
effect essentially does not reject any event detectable by 
GOES. 

Microflares analyzed in a previous RHESSI spectral 
study of flares of all scales, 42 microfl ares out of a sam- 
ple of 85 flares, (|Battaglia et al.l 120051 ) are shown as the 
crosses in Figure [13] and also showed no correlation. 
These events specially chosen to cover a wide range of 
RHESSI Q^axe magnitudes were analyzed during the peak 
in 12-25 keV. This is earlier in the flare phase than for our 
survey. Thus they have correspondingly higher tempera- 
tures and lower emission measures. Anothe r recent study 
of 18 microflares from a single active region fStoiser et alj 
l2007f) found similar results to the Battagiia ct al. (2005)) 
study. 

The dashed line in Figure [13] is the correlation 
found using soft X-ray observatio ns with GOES and 
Yohkoh/BGS (iFeldman et al.lll996bD . for aU size of flares 
from A through to X class, not just microflares. This 
correlation has a order of magnitude spread in emis- 
sion measure and only becomes apparent when a large 
range of flare magnitudes are included; for only small A- 



FlG. 16.— Rff£;55/ non-thermal photon flux at 12 keV, from the 
fitted broken power-law, against fitted temperature (top) and emis- 
sion measure (bottom). The correlation plots are the microflares in 
this study, the flt to this data is the dot-dashed line. The dashed 
line is fro m a previous RHESSI study using the non-thermal flux 
at 35 keV l|Battaglia et al.ll2005l V 

class events the correla tion was in the opposite direction 
(jFeldman et alJ Il996al) with the emission measure de- 
creasing with increasing temperature. In comparison the 
i?ffi?5'iS'/ data, which i s observed at higher energies than 
iFeldman et al.l ()1996bD . show higher temperature and/or 
lower emission measures. This is consistent with the 
DEM peaking at temperatures lower than those observed 
with RHESSI and closer to those lower temperatures ob- 
served by GOES and Yohkoh/BGS. This was found to be 
case for the DEM of a large flare observed in soft X-rays 
with GOES an d Yohkoh/SXT and in har d X-ra ys with 
Yohkoh/HXT (jAschwanden fc Alexande^ l2001h . For 
RHESSI emission measures above 8 x lO'*^ cm~'^ there is 
the hint of a similar correlation with temperature but it 
is obscured by the temperature selection effect at lower 
emission measures. Further studies of solar flare temper- 
ature and emission measure suggested that the emission 
measure of EUV nanofla res approximately scales as 
(jAschwanden et al.|[2000l ). but using various studies over 
larger scales suggested th e emission measure may scale 
as T"' (jAschwandenI l2007l ) . These scalings are only ap- 
proximate as there is large scatter about the correlation 
line. Looking over larger ranges with different instru- 
ments helps to reduce the overall influence of the selec- 
tion effects but there is still ambiguity as to how these 
parameters scale. 
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Fig. 17.— RHESSI non-thermal flux at 12 keV, from the fitted 
broken power law, against the the thermal flux over 4-8 keV, from 
the thermal model. The dot-dashed line is the fltted correlation 
line. 

To gain a better understanding of the RHESSI tem- 
perature and emission measure, we have also calculated 
the GOES temperatures and emission measures for each 
of these microfiares. This was done using the same pre- 
flare time for background and peak time of emission in 
6-12 keV as was used in the RHESSI analysis. This was 
possible for 6,740 microfiares, removing those events for 
which a GOES temperature and emission measure were 
not reliably calculable. The resulting correlation plots 
are shown in Figure [Ml The bottom right hand of this 
correlation plot (the lowest GOES temperature for each 
RHESSI temperatures) does seem to be directly propor- 
tional, except the i?ffi?5'S'/ temperatures are about 5 MK 
higher. The rest of the plot is again dominated by the 
selection effect in the i?ffi?5S'/ temperatures. The events 
with GOii'S' temperature < 10 MK have been highlighted 
by the white contour in Figure 1131 If the assumption is 
that the RHESSI temperate estimate is too high in these 
events, and the emission measures is consequently too 
small, then this set of highlighted events should be moved 
upwards and to the left in Figure [T31 This shift produces 
a clearer hint of the previously found correlation between 
temperature and emission measure. 

The emission measures in Figure [14] do not show any 
such temperature selection effect and are nearly directly 
proportional, with the GOES emission about twice that 
observed in RHESSI Again this will be due to RHESSI 
observing in a temperature range which is higher than 
the peak temperature in the DEM. The proportional- 
ity between the RHESSI and GOES emission measures 
on face value suggests a similarly steep DEM in all the 
microfiares but this might be arising from the relative 
instrumental sensitives. It is important to remove these 
instrumental effects to recover the underlying DEM but 
this is a complicated process a nd has only been successful 
for in dividual large flares, (c.f. lAschwanden fc Alexander! 
[200l . ^ 

Ot her studies (jFeldman et al.l Il996bl : iBattaglia et al.l 
12005^) have also investigated how the temperature and 
emission measure relate to the GOES flux of the 
events. The RHESSI microflare temperature and emis- 



sion measure plotted against each event's correspond- 
ing background-subtracted GOES l-SA is shown in Fig- 
ure [ini Again there is a hint of a correlation between 
the RHESSI temperature and GOES flux for the events 
above B-class (10^^ Wm~^), but it is obscured in smaller 
events again due to the temperatu re selection effect. The 
previous study of RHESSI &ares (|Battagha et al.ll2005f) 
found a correlation when using large flares in addition 
to microfiares. This correlation and the microfiares in 
their sample are shown in Figure [T51 by the dashed line 
and crosses, a nd is steeper than fou nd using GOES and 
Yohkoh/BGS (jFeldman et al.lll996bD . The hint of a cor- 
rel ation in our microflare study s cales in a manne r closer 
to iFeldman etall (|1996b[ ) than IBattaglia et"all 12005), 
although the temperatures found are consistently higher. 

There is a clear correlation between the emission mea- 
sure and GOES flux, which can be fitted as Fq — 
1.15xlO-52£;^/0.96^ ^^^^ Wm-2 a nd EM in cm'^ 

A sim ilar result was also found in the IBattaglia et al.l 
(1200^ study, fitted over a large flare to microflare range 
finding Fq = 3.6 x 10-'='°EM°-^'^ . The consistently lower 
emission measures in this study are again due to this 
analysis occurring earlier in the microflare than in our 
study. 

We can also investigate how the non-thermal emis- 
sion relates to the thermal parameters, which is shown 
in Figure 1161 Here the non-thermal flux i^N , in units 
of photon flux at Earth s^^cm^'^keV^^, is taken from 
the fitted broken power-law parameters at 12 keV. This 
energy is used as it is above the break energy ee in 
all events but is still at an energy for which we ob- 
serve non-thermal emission. Both the temperature and 
emission measure scale with the non-thermal emission. 
The correlation with the temperature is fitted as T = 
1.41 log Fn -I- 13. 2, with T in MK. A si milar correlation 
was found in the IBattaglia et all ()2005[ ) study using the 
flux at 35 keV for the non-thermal flux, since larger flares 
were also analyzed, finding T = 1.46 log ^35 -I- 21.57. We 
find that the scaling of the emission measure is flatter, 
EM = 4 X io46f 0.57^ compared to t he previous study 
EM = 5 X lO^SFgOgSi (|Battagha et al.l[2005h . again EM 
in units of cm~'^. 

In FigurefTTIthis non-thermal photon flux at 12 keV F-^ 
is plotted against the model thermal flux over 4-8 keV, 
Ft. The thermal and non-thermal emission correlate, 
with the fit being Fy, = 3.4 x lO'^F^^ ^^, but with a 
greater scatter in the smaller events. At higher fluxes 
the correlation steepens, suggesting that there might a 
greater proportion of non-thermal emission relative to 
thermal emission in larger flares. However, this may just 
be an instrumental effect, arising from detector pileup 
and livetime issues in these larger events prior to the 
shutters deploying. 

4. ENERGY DISTRIBUTIONS 

4.1. Thermal Energy Frequency Distributions 

Using the volumes found in Sj2| and temperature and 
emission measure found in ^J3| the thermal energy Wt , 
over the time of peak emission in 6-12 keV, can be calcu- 
lated via equation ([3]). The energies range from 10^^ erg 
to 10"^" erg with the median being about 10^® erg. The 
frequency distribution (number of events per energy bin 
range, area of solar disk and duration of observation pe- 
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Fig. 18. — Frequency distribution of the thermal energy, of 9,161 RHESSI microfiares (RH) in context of thermal energy distributions of 
nanoflares (TA, TP and EB) and active region transient brightenings (SS). The dotted line indicated a power-law index of a = 2. 
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Fig. 19. — Frequency distribution of the microflare thermal en- 
ergy, each line indicating a different year. The number quoted 
next to the year is the number of events from each year used in the 
distribution. 



riod) of this energy for 9,161 microfiares is shown in 
Figure [TH] The resulting RHESSI thermal energy dis- 
tribution is not a clear power-law; it has a turn-over at 
low energies and steepens at higher energies. These fea- 
tures are instrumental effects due to missing the smallest 
events, with insufficient counts to either find the events 
or successfully analyze them, and the largest events, due 
to detector livetime issues before the attenuating shut- 
ters come in. Considerably more small events are miss- 
ing than large, since the discrepancy from a power-law 
is greater at lower energies than high. As this distribu- 
tion deviates from a power-law it will not be fitted to 
obtain the power-law index a. Instead, Figure [T51 shows 
an a = 2 line, indicating the parts of the RHESSI ther- 
mal distribution can be steeper and flatter, or larger and 
smaller, than a = 2. 
Compared to the previous distributions found for EUV 



nano f lares (iKrucker fc Bend 119981: lAschwanden et al.l 
[2OOOI: iParnell fc Juod 120001: iBenz fc Krucked I2002D and 
soft X -ray active-region transient brightenings (jShimizul 
mil), the RHESSI energy distribution appears as an 
extension to these at higher energies. This is both re- 
markable and deceptive since these distributions were 
found for very different types of events, using various 
instruments and for different periods d uring the solar 
cycle. For instance the SXT energies jSh imiz 
are from 291 brightenings in one active re gion over 5 
days in August 1992, whereas the EIT (.Krucker fc Bend 
ll998HBenzlrKr ucker 20021 and TRACE'EVY quiet sun 
obser vation ([Aschwanden et al.l 120001 : IParnell fc JuppI 
[2OOOI) were found over about an hour each on 12 July 
1996, 17 February 1998 and 16 June 1998 respectively. 
So there are two key issues that have to be taken into 
account when looking at the energy distributions in Fig- 
ure 1181 First, various instruments were used, so differ- 
ent components of the thermal energy will be observed 
and with distinctive instrumental selection effects will 
influence each distribution. This makes it difficult to de- 
termine whether these are similar events or completely 
distinctive physical processes. Second, the distributions 
cover different phases of the solar cycle. The previous 
studies show a snapshot of the energy distribution of dif- 
ferent small energy release features in the solar corona 
whereas the RHESSI microflare energy distribution rep- 
resents 5 years of observations of the declining phase of 
the solar cycle. The difficulty here lies in determining 
whether these snapshot surveys demonstrate typical or 
unusual behavior. Using the RHESSI thermal energies 
we can investigate how the energy distribution changes 
with time by plotting the distributions for each year 
separately, as shown in Figure [TH These distributions 
have similar shapes except that the normalization de- 
creases by over an order of magnitude between 2002 and 
2006. So the non- RHESSI distributions in Figure [H] 
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could be shifted vertically by a considerable amount if 
they were found during a different part of the solar cy- 
cle. This invariance in the shape of the distribution dur- 
ing t he solar cycle has been found previous ly in both 
soft (Fcldman et al.' 1997; Vcronig ct an i2002l ) and hard 
(fCrosb y et al.„ .1993.: ,Lu et al.. 1993) X-rays, and also il- 
lustrates that the selection effects on the RHESSI data 
do not vary with time. 

4.2. Non-Thermal Power Frequency Distributions 

To calculate the power in accelerated electrons from 
the information about the power-law in the photon spec- 
trum, we use 7 and Jq (the power-law index and normal- 
ization), found in Sj3l in equation ([4|). However there is 
ambiguity as to the low-energy cut-off Eq because the 
observed photon spectrum depends only weakly on it. 
In addition the RHESSI imcmQ.aie spectrum covers both 
the thermal and non-thermal energy ranges and so there 
is inherent ambiguity. Previous studies used their X-ray 
threshold for E g resulting in large uncertainties in th e en- 
ergy estimates (jCrosbv et al.lll993HLrn et al.ll2001f ). As 
RHESSI makes spectral measurements down to the ther- 
mal component we can provide a better estimate of the 
non-thermal energy content. In Part I of these papers 
a fixed Eq was used to provide a r ough estimate of th e 
non-thermal power at peak time (jChriste et al.l [20081 ). 
The full spectrum fitting in ^provides a better estimate 
of the parameters required to calculate the non-thermal 
power via equation (jj]). However, despite this improve- 
ment we only have an estimate of where the photon spec- 
trum begins to flatten (ee) and not the actual cut-off in 
the electron distribution Eq. 

To obtain an estimate of Eq we have therefore investi- 
gated how the fitted broken power-law model relates to 
the expected photon spectrum from the electron distri- 
bution using the numerical integration code of Holmam 
(|20G3). An example of this is shown in the left panel 
of Figure 1201 where the photon spectrum for an electron 
distribution with (5 = 8 and Eq was calculated and then 
fitted with a broken power-law in the same manner as was 
used for the microflare spectrum. T his broken po wer-law 
has 7 = 7.03, expected as i5 = 7 -I- 1 (|Brownlll971h . above 
a break of ee = 6.06 keV. Repeating this process for var- 
ious values of the electron distribution parameters 5 and 
Eq^ shows how the resulting fitted e-Q scales with these 
parameters. This is shown in the middle panel in Fig- 
ure [501 where ee is given against Eq for three different 
values of 5. For a single value of 5, Eq approximately 
scales linearly with ee and this scaling steepens with in- 
creasing 5. So we can approximate this relationship to 
first order by linearly fitting the relationship between eb 
and Eq for various 5 and then linearly fitting how the pa- 
rameters of the first fit vary with 6. This empirical rela- 
tionship between the observed parameters of the photon 
power-law 7, eb and the low energy cut-off of the electron 
distribution Eq can be found: 

Eq k 0.157 -f (1.86 - 0.047)eB - 3.39. (6) 

Using equation ([6]) the values of 7 and eb for the 4,236 
microflares with trustworthy non-thermal spectral fits re- 
sult in a histogram of Eq for these events, shown in the 
right panel in Figure 1201 The low energy cut-offs range 
from 9 to 16 keV with the median being about 12 keV, 
with generally eb ~ 0.75i?c- The uncertainty in Eq can 



be seen in the average ratio of the error (found from the 
fit errors in 7 and eb) to Eq also shown in the right 
panel of Figure [201 For Eq < 12 keV, where the major 
of the event lie, the uncertainty is about 10% and drops 
to about 1% for the few events with larger Eq. 

These values of Eq can then be used in equation ([4]) 
to calculate the power in non-thermal electrons above 
Eq, P(> Eq), the frequency distribution of which is 
shown in the left panel in Figure [211 Here the power 
ranges over 10^^ to 10^^ erg s~^ with the median be- 
ing about 10^^ erg s~^. Also shown is the ratio of 
the error to the power, as a function of the power. 
The few microflares which show a power of around 
10^® erg s~^ have the smallest errors with uncertainties 
about 10%. These powers seem relatively high for small 
flares, as the RHESSI power es timat es in large flares are 
10^^ to lO'^" erg s" ^ ( Holm a^ [2001 Emslic et'al][200l 
ISaint-Hilaire fc Benz>,2005 : Sui et al.ir2005, . ,200:i7though 
it should be noted that the non-thermal emission in mi- 
croflares lasts for « 10 seconds whereas it can last for tens 
of minutes in large flares. The total non-thermal energy 
content in large flares is many orders of magnitude larger 
than in microflares. The majority of the microflares 
shown in Figure [^H show non-thermal powers consider- 
ably smaller than this level, although with increasing un- 
certainty. The median power of 10^^ erg s~^ has about 
a 50% error and the smallest events at 10^^ erg s~^ have 
almost 100% error. It is this large uncertainty in the 
power that results in this distribution deviating from a 
power-law. Although the selection effects and biases will 
affect this non-thermal distribution as seen in the ther- 
mal distribution, these effects are hidden by the large 
uncertainties. Since it deviates from a power-law, it has 
not been fitted using this model. 

To allow comparison to the power distribution in elec- 
trons above 25 keV P(> 25) for 2,878 large flares fo_und 
over 1980 to 1982 using SMM/HXRBS dCrosbv et all 
[199I the same power is estimated in the RHESSI mi- 
croflares, shown in the right panel in Figure [21] These 
powers using Eq = 25 keV for the RHESSI microflare 
ranges over lO'^^ to 10^^ erg s^^ with the median power 
being about 10^^ erg s~^. On this basis the RHESSI mi- 
croflares cover events down to two orders of magnitude 
smaller than those in the SMM/HXRBS study. Again 
the error is smallest in the largest power estimates. This 
P(> 25) i?_ff£'S'5'/ distribution deviates less from a power- 
law than the P{> Eq) distribution. This is because there 
is less uncertainty in the power estimate in the larger 
events as the fixed Eq = 25 keV has no associated error. 

5. DISCUSSION & CONCLUSIONS 

We have analyzed 25,705 microflares, successfully re- 
covering trustworthy spatial information about the ther- 
mal emission for 18,656 of them, as well as the thermal 
spectral fit in 9,161 and the thermal and non-thermal fit 
in 4,236. The median values and ranges of each of the 
microflare para meters are summ arized in Table [1] As 
found in Part I (jChriste et al.ll2008.) the microflares occur 
only in active regions and this immediately suggests that 
these events are ill-suited to heat the overall corona. It 
also helps to show that the X-ray emission outside active 
regions, especially during periods of quiet Sun, is consid- 
erably smaller than the smallest active-region flares ob- 
served with RHESSI. This is conflrmed by the RHESSI 
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this would be complicated by systematic effects related 
to RHESSFs attenuating shutters. An alternative would 
be to investigate this selection effect by comparing model 
spectra for known temperatures and emission measures 
with those found by fitting this spectrum once noise 
and the instrumental response have been included. This 
has alr eady been attempted for the TRACE nanoflares 
study (lAschwanden fc Parnel]|[2002f l. and it revealed a 
difference between the observed parameters and the "in- 
trinsic" ones. Another way of investigating the instru- 
mental response and biases is in comparison to results 
from other instruments. We have found that the GOES 
emission measures are typically a factor of 2 larger than 
those found with RHESSI. This maybe due to the un- 
derlying DEM in microflares scaling similarly but might 
also be due to the instrument's relative sensitivity. This 
could be further studied by forward-fit modeling of the 
instrumental r esponses to recover the underly ing DEM 
in each event (jAschwanden fc Alexander! 1200 If ), instead 
of assuming a constant emission measure. 
The thermal energies for the peak time in these mi- 



quiet Sun study (|Hannah et al.|[2007f) which found a limit 
on the fiux over two orders of magnitude smaller than 
that of the smallest microfiares found here. 

Forwarding-fitting the visibilities of RHESSI data is a 
fast and efficient wa y to recover the spatial information 
(jHurford et al.ll2005| ). We find that the thermal emission 
(4-8 keV) of the microflares shows predominantly loop- 
like structures that have a median width of 11" and 32" 
length. Those events where no spatial information was 
recoverable were predominantly those with the fewest 
counts, although some larger events produce poor spa- 
tial information as there were other instrumental issues. 
These spatial scales do not correlate with the magnitude 
of the microflare, the 4-8 keV flux from the loop nor 
the GOES class, and so small flares are not necessarily 
spatially small. This may be due to the prevalence of a 
typical loop scale size associated with active regions, and 
it is independent of the amount of energy the fiare has 
deposited into heating the material that evaporates to fill 
these loops. In the largest fiares this energy might be de- 
posited over a larger area, evaporating material into more 
loops, hence the larger structures and arcades that are 
observed. The widths of these loops may not be as well 
resolved as the lengths, with the observed single RHESSI 
loop possibly being several narrow and long loops be- 
side each other, as observed in microfi are with RHESSI 
and Hinode/XRT (jHannah et al.ll2008f l. The volume of 
these loops can be estimated, the median value being 
1 X 10^^ cm'^. This volume estimate assumes a cylindri- 
cal geometry and filling factor of unity, meaning that our 
volumes, and hence thermal e nergies, are upper limits, as 
this filling factor can be^l (|Cajgill fc Klimchuk|[l997l : 
iTakahashi fc Watanabdl2000f ). 

We find that microflares are hot, with a median tem- 
perature of 13 MK. In these microflares the median emis- 
sion measure is 3 x 10*^ cm~'^, which combined with 
the volume allows a density of 6 x 10^ cm~'^ to be cal- 
culated. As the volume may be an overestimate, this 
means that the density might be underestimated. The 
presence of the spectral feature due to the Fe K-shell 
transitions proves that there rea lly is hot plasma present 
with T > 8 MK (|Phillipsll200l . Correlations between 
the temperature and other parameters are difficult to ex- 
tract from our sample due to the temperature selection 
effect, resulting in only a relative narrow range of tem- 
peratures being sampled. Studying larger fiares would 
aid by extending the range of flare temperatures, but 



crofiares were found to range from 10 ° erg to 10 erg 
with the median energy being 10^* erg. The thermal en- 
ergy distribution deviates from a power-law at low and 
high energies but this can be explained by selection ef- 
fects and does not suggest that the underlying true distri- 
bution is not power-law. The smallest events are missing 
as they do not have enough counts to be found either 
above solar or instrumental background, successfully im- 
aged or spectrally analyzed. The larger events are miss- 
ing from the distribution due to RHESSFs attenuating 
shutters, so either they are exclude from our selection 
or have poor detector livetime. Although it is possible 
to extend this energy distribution to higher energies by 
analyzing all large fiares observed with RHESSI, there 
will still be a strong instrumental effect at the shutter 
transition due to detector livetime effects. The statisti- 
cal errors in the fitted parameters will affect the thermal 
energy estimates but these have a relatively small effect 
compared to the selection effects. The greatest uncer- 
tainty is the systematic error that arises from the filling 
factor used to estimate the therma l volume: taking an 
extre me filling factor of f ~ 10~^ (jCargill fc KlimchukI 
Il997t l would have the effect of reducing the energies by a 
factor of 100. Therefore the thermal energies quoted here 
about the time of peak emission in 6-12 keV are upper 
limits. The comparison of the RHESSI thermal distribu- 
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tion to other the thermal distributions of other transient 
coronal energy releases is difficult as these events are from 
a limited sample of events over a very short time period. 

With RHESSI we are able to investigate the flattening 
of the non-thermal photon spectrum and can empirically 
estimate the low energy cut-off Ec in the electron distri- 
bution from the observed photon spectrum. However we 
are only able to successfully fit the non-thermal spectral 
component, in addition to the thermal component, in 
4,236 events. The difficultly in successfully finding and 
fitting this non-thermal emission arises from the steep 
spectra that start close to the thermal spectral features 
about 7 keV. So, although the non-thermal component 
was only successfully fitted in a minority of events, it does 
not mean that there is a lack of accelerated electrons in 
the others. It does mean however that the relatively large 
uncertainties in the non-thermal parameters will result 
in large uncertainties in Ec and the non-thermal power 
estimates. These large uncertainties and the biases intro- 
duced by rejecting many events due to poor spectra fits 
can be easily seen in the histograms of the non-thermal 
parameters 7 and cb- We find that the microfiares typi- 
cally have non-thermal emission represented by a broken 
power-law with index 7 = 7, break at 9 keV. We can es- 
timate Eq — 12 keV, for breaks in the photon spectrum 
above 7 keV. 

These non-thermal parameters were used to calcu- 
late the power in accelerated electrons above Ec, which 



range over 10 to 10 erg s~ with a median value of 
10^^ erg s^^. This distribution again deviates from a 
power-law at low and high energies due to the uncertain- 
ties in the power estimates and the rejected events due to 
poor spectra fits rather than the selection effects present 
in the thermal distribution. The uncertainties in the few 
events with the largest non-thermal component are only 
~ 10% and so this cannot explain the fact that in some 
of these small flares the rate of energy rel ease is compa- 
rable to larger flar e s. 10^"^ to 10^° erg s~^ (lHolinanll200l 



Emslie et al.l l2004t ISaint-Hilaire fc Bendl2005l : ISui et al l 
2005, 2G07|. However the non-thermal emission in mi- 
croflares lasts for only w 10s seconds whereas it can last 
for tens of minutes in large flares. The non-thermal en- 
ergy content in large flares is thus many orders of mag- 
nitude larger than in microflares. To make a direct com- 



parison of these RHESSI microflare results to the peak 
non-thermal power distribution found for large flares by 
ICrosbv et all ([1993), the same instrumental Eq — 25 keV 
has to be used. On this basis the i?_ffi?5'S'/ microflare peak 
pow ers extend to two or ders of magnitude smaller than 
the (ICrosbv et all I1993D s tudv. The true non-thermal 
power in the ICrosbv et ah! ()1993f ) flares cannot be deter- 
mined as the spectrum was not observed to sufficiently 
low energies and so that study could have easily under- 
estimated the non-thermal energy content. We conclude 
that the instantaneous non-thermal power in a microflare 
can be surprisingly large. 

One further explanation for the unexpectedly large 
non-thermal power in the microflares might be that the 
physical model used is less suited for small flares than 
large. The energy in the accelerated electr ons was foun d 
using the standard cold thick target model (|Brownll971h . 
It has been suggested that for the lowest energy elec- 
trons the target is warm and not cold, and 5kT could 
be used as an approximate cut-off ( Emslie ,2003) . Un- 
fortunately for the microflares, this gives an even lower 
cut-off at around 4-5 keV, which results in a even larger 
non-thermal power estimate. Nevertheless the idea of the 
non-thermal electron distribution smoothly transitioning 
into the thermal distribution seems physically more re- 
alistic than a sharp low energy cut-off and may provide 
a clue to the physics of the energization of the electron 
distribution function. Another change to the model of 
the hard X-ray ei nission wou ld be the inclusion of free- 
bound emissions (jBrown fc M allik 200J,), as our models 
currently only assume free-free continuum. Microflares 
might be the ideal type of flares to study this seldom- 
considered mechanism as the resulting additional spec- 
tral features would occur above 10 keV and would likely 
have been hidden by the thermal component in large 
flares. However the temperatures requi red for such free- 
boun d features to be present, > 20 MK (jBrown fc MallikI 
l2007f) , are not typically observed in RHESSI microflares. 

The analysis presented here is only performed at the 
time of peak 6-12 keV emission even though the most 
desirable results would investigate the full time range of 
each microflare, producing the total thermal and non- 
thermal energy estimates. But as this paper has shown, 
it is a considerable undertaking just analyzing this peak 
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TABLE 1 

Median and range of RHESSI microflare parameters. 



Microflare Parameter 




Median Value 


Range 


(5% to 95%) 


Duration^ 


D 


5.4 mins 


2.2-15 mins 


Temperature 


T 


12.6 MK 


10.7 


- 15.5 MK 


Emission Measure 


EM 


3 X 10-*^ cm-3 


4 X 10-*5 - 


- 2 X 10'*'' cm-3 


Thermal Loop Width 


w 


8 Mm (11") 


3-20 Mm (4 - 28") 


Thermal Loop Length 


I 


23 Mm (32") 


7-77 Mm (10 - 107") 


Thermal Loop Volume 


V 


1 X 10^'^ cm^* 


5 X 10^5 


- 2 X 1028 cm3 


Density 


ric 


6 X 10^ cm"^ 


8 X 10* - 


- 3 X 10^° cm-3 


Power-law Index 


7 


7 




4-10 


Break Energy 




9 keV 


7 


- 12 keV 


Low Energy Cut-off 


Ec 


12 keV 


9 


- 16 keV 


Thermal Energy'' 


Wt 


10^** erg 


1026 


- 10^° erg 


Non-thermal Power'' 


Pn(> Ec) 


10^6 erg s-l 


1025 _ 


10^8 erg s~' 


Non-thermal Power' 


Pn{> 25) 


10^* erg s~' 


1022 - 


lO^'^ erg 



Values from part I of this article dChriste et al.|[2008l ): 1 min lower limit due to selection 
effects.' All parameters are estimated from analysis using 16 seconds around the time of 
peak emission in 6-12 keV, so resulting energy estimates for those times. 



time period. The next steps in this microflare study 
would be to investigate the biases and instrumental ef- 
fects using simulations and making improvements to the 
forward-fit model of the emission. Although RHESSIha.s 
for the first time allowed the thermal energy and non- 
thermal power distributions to be studied systematically 
in microflares, the uncertainty in the transition between 
thermal and non-thermal components highlights the in- 
strumental effects and biases. We also have suggested 
the possibility that the standard hard X-ray model used 
is wrong. The model and instrumental effects could be 
studied by investigating the microflares not only with 
RHESSI but other instruments as well. It may however, 



require an instrument with better energy resolution and 
lower instrumental background than RHESSI before we 
can unambiguously determine the energy component of 
microflares and decide whether there is an issue with the 
hard X-ray model used. 
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